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Abstract 

Recently different formulations of the first-order Lighthill-Whitham-Richards (LWR) 
model have been identified in different coordinates and state variables. However, there 
exists no systematic method to convert higher-order continuum models into car-following 
models and vice versa. In this study we propose a simple method to enable systematic 
conversions between higher-order continuum and car-following models in two steps: 
equivalent transformations of variables between Eulerian and Lagrangian coordinates, 
and finite difference approximations of first-order derivatives in Lagrangian coordinates. 

With the method, we derive higher-order continuum models from a number of well- 
known car-following models. We also derive car-following models from higher-order 
continuum models. For general second-order models, we demonstrate that the car¬ 
following and continuum formulations share the same fundamental diagram, but the 
string stability condition of a car-following model is different from the linear stability 
condition of a continuum model. This study reveals relationships between many existing 
models and also leads to a number of new models. 

Keywords: Eulerian coordinates; Lagrangian coordinates; Primary models; Secondary 
models; First-order models; Higher-order models; Continuum models; Car-following models; 
Fundamental diagram; Stability. 


1 Introduction 


In the three-dimensional representation of traffic flow (Makigami et ah, 1971), as shown in 
Figure [lj the evolution of a traffic stream can be captured by a surface in a three-dimensional 
space of time t, location x, whose positive direction is the same as the traffic direction, and 
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Figure 1: The x — t — N surface for a traffic stream 


vehicle number N, which is the cumulative number of vehicles passing location x at t after a 
reference vehicle (Moskowitz, f965). We call these primary variables of traffic flow. In the 
figure, the red curve, 


N = n(t, x ), 


(la) 


represents the x — t trajectory of vehicle N, which can be collected with GPS devices; the 
green curve, 


x = X{t,N), (lb) 

represents the cumulative number of vehicles N before t at location x, which can be extracted 
from vehicle counts of a loop detector; and the blue curve, 

t = T(x,N), (lc) 

represents the number of vehicles N downstream to x at time t, which can be counted from 
a camera snapshot. Among the three primary variables, the two independent variables form 
a coordinate system and the dependent one is the state variable. In particular, (£, x) form 
the Eulerian coordinates, and (t, N ) the Lagrangian coordinates. 

The goal of a traffic flow model is to determine the surface from given initial, boundary, or 
internal conditions. It is usually written as a partial differential equation or a corresponding 
difference form. We call a traffic flow model primary if one of the three primary variables 
is the unknown state variable, and the order of the model is the maximum order of partial 
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derivatives. In contrast, if a partial derivative of a primary variable is used as the unknown 
state variable, we call the model secondary. In this study, we will categorize traffic flow 
models according to their coordinates, state variables, and orders. 

In the Eulerian coordinates, we can define traffic density k = —n x and flow-rate q — n t , 


where the subscripts mean partial derivatives. Then the celebrated LWR model (Lighthill 


and Whitham, 1955; Richards, 1956) is a partial differential equation in k: 

= 0, 


dk d(j)(k ) 

dt 


dx 


( 2 ) 


which is a hyperbolic conservation law (Lax, 1972). Thus, the LWR model is a first-order 


secondary model in the Eulerian coordinates, abbreviated as the E-S-l model. The model 
can be derived from the following three rules: 

• Fundamental law of traffic flow: speed v — |: 

• Conservation of traffic flow: kt + q x = 0; 

• Fundamental diagram of flow-density relation: 

q = <f>(k). (3) 

In particular, the following triangular fundamental diagram has been shown to match 
observations in steady traffic well (Munjal et al. 1971 |Haberman 1977 Newell, 1993): 


c t>(k ) = mm{vfk,w(kj — k)}, 

where Vf is the free-flow speed, — w the shock wave speed in congested traffic, and kj the 
jam density. The Greenshiclds fundamental diagram is (Greenshiclds, 1935) 


(f)[k) = Vfk( 1 — k/kj). 


In (Newell, 1993), the E-P-l model, i.e., the first-order primary traffic flow model in the 


Eulerian coordinates, was first discussed: 


n t - cj)(—n x ) = 0, 


(4) 


which is a Hamilton-Jacobi equation (Evans 


1998). This model can be numerically solved 


with the variational principle (Daganzo, 2005a[|b), or the Hopf-Lax formula (Claudel and 


Baycn, 

2010 

£ 

HH 

Daganzo, 

2006b) 


first-order primary formulation in the Lagrangian coordinates, was derived based on the 
inverse function theorem as 


X t .~0{-X N ) = 0, 


(5) 
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Figure 2: An observed flow-occupancy relation (Hall et al. 1986) 


which is also a Hamilton-Jacobi equation. If we denote the spacing by S' = —Xn, then the 
speed-spacing relation is 


0(S) = S<p(l/S). 


( 6 ) 


In particular, for the triangular fundamental diagram 

9(S) = mm{vf,w(kjS — 1)}; 
and for the Greenshiclds fundamental diagram 


6(S) = v f {\ 


Sk 


-)■ 


Further in (Daganzoj, 2006b ja), various discrete forms of (j5j) was obtained based on the 
variational principle. In particular, Newell’s car-following model is shown to be a discrete 
form of ([5]) with the triangular fundamental diagram. In (Leclercq et al., 2007), the L-S- 
1 formulation was derived as a hyperbolic conservation law and numerically solved with 


the Godunov method (Godunov, 1959). More recently, in (Laval and Leclercq, 2013), the 


equivalent primary and secondary formulations in the (x, IV)-coordinates, called T-models, 
were derived and solved. In summary, all the primary formulations of the LWR model are 
Hamilton-Jacobi equations and can be solved with the variational principle, the Hopf-Lax 


formula, or other techniques (Evans, 1998); and all the secondary formulations are hyperbolic 


conservation laws and can be solved with the Godunov scheme. 

An equilibrium flow-density relationship (|3| has been derived from various car-following 


models in steady states (Gazis et al. |1961 Haberman 1977 Chapter 61) and verified with 


empirical observations (Del Castillo and Benitez 1995 Cassidy, 1998[). However, in reality, 
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traffic can be in non-equilibrium states, and the flow-density or flow-occupancy relations are 
not unique, as shown in Figure [2} In such non-equilibrium traffic, instability and hysteresis 


phenomena have been observed (Newell, 1965 Treiterer and Myers, 1974 Sugiyama et al. 


2008). In order to capture such non-equilibrium features of traffic flow, researchers have 


introduced numerous higher-order continuum models: in (e.g. 

Payne, 

1971 

Whit ham 

1974 


Zhang, 1998; Aw and Rascle, 2000; Zhang, 2002), the conservation equation is coupled with 
one or more equations in speed or other variables. There have been extensive discussions 


on theoretical consistency of some higher-order models ( 

Daganzo, 

1995 

Lebacque et al. 

to 

o 

o 

-u 

Helbing and Johansson, 

2009; 

Zhang, 

2009) 

Helbing 

5 , 2009) and on the stability and 

other properties ( 

Kerner and Konhauser| 

1994 

Li, 

2001; 

fin and Zhang, 

2003 

). Note that 


however, these models are still limited, as they are just capable of explaining a subset of 
many non-equilibrium phenomena (Kerner, 2004). Non-equilibrium traffic phenomena have 
also been successfully modeled with microscopic car-following models since 1950’s in the 
transportation field ( |Pipes 1953; Herman et al. 1959) and in the physics field, as traffic 
systems can be considered many-particle systems (Bando et al. 1995). Recently stability and 


other properties of car-following models have been subjects of some mathematical studies 


(Wilson 

2001 

Orosz and Stepan, 

2006 

). Reviews of various car-following models are available 

in ( 

Rothery 

1992; 

Brackstone and McDonald, 

1999; 

Chowdhury et al. 

2000 

Helbing 

2001) 


In the literature, many higher-order continuum models have been derived from car-following 


models using the method of Taylor series 

expansion (e.g. I 

3 ayne, 

1994; 

Zhang 

1998; 

Berg et ah, 

2000; 

Zhang 

2002; 

Jiang et al. 

,2002) 


1971: Del Castillo et al. 


process is model specific, and different approximation methods can lead to different continuum 
models for the same car-following model. In addition, such a method cannot be used to derive 
car-following models from continuum models. 

2006a,b 


From the studies on different formulations of the LWR model in (Daganzo 


Leclercq et ah, 2007 Laval and Leclercq, 2013), we can see that car-following models can be 


considered as discrete primary formulations in the Lagrangian coordinates, and continuum 
models as secondary formulations in the Eulerian coordinates. Based on this observation, 
in this study we will develop a systematic method for conversions between higher-order 
continuum models and car-following models. However, the framework is different from 
that for the first-order models: first, the higher-order L-P models are no longer Hamilton- 
Jacobi equations, and the variational principle or the Hopf-Lax formula cannot be applied to 
obtain car-following models from the L-P models; second, as higher-order derivatives of the 
primary variables are involved, the inverse function theorem cannot be used to transform 
variables between Eulerian and Lagrangian coordinates. To resolve the first difficulty, we 
use simple finite difference methods to derive car-following models from L-P models, since a 
driver’s car-following behavior is only impacted by its leader; to resolve the second difficulty, 
we systematically derive the transformations of related variables between Eulerian and 
Lagrangian coordinates. The conversion method consists of two components: finite difference 
approximations of L-P derivatives, and transformations between E-S and L-P variables. 
With this method, we will derive continuum formulations of car-following models and vice 


5 






































































































































































































versa. We will be able to reveal relationships between some high-order models and also 
derive new continuum and car-following models from existing ones. We will also examine 
the equivalence between general second-order continuum and car-following models in their 
steady-state speed-density relations and stability of steady states. 

The rest of the paper is organized as follows. In Section 2, we present a method for 
conversions between continuum and car-following models and apply it to the first-order 
models. In Section 3, we derive continuum formulations of some car-following models. In 
Section 4, we derive car-following formulations of higher-order continuum models. In Section 
5, we examine the equivalence between continuum and car-following models. In Section 6, 
we discuss some future directions. 


2 A systematic conversion method 

In this section we present a systematic method for conversions between car-following and 
continuum models. First, we present transformations between E-S and L-P variables, which 
lead to equivalent E-S and L-P formulations of a traffic flow model. Second, we present finite 
difference approximations of L-P. derivatives 


2.1 Transformations between E-S and L-P variables 


In E-S models, variables include k , q, v, and their derivatives with respect to x and t ; in L-P 
models, variables include X and its derivatives with respect to N and t. In this subsection, 
we establish the equivalence between E-S and L-P models through transformations between 
Eulerian and Lagrangian variables. 

For a physical variable, we denote it by F(t,N ) in the Lagrangian coordinates and by 
f(t,x ) in the Eulerian coordinates. Then from (fla|) we have 



f(t,x) = F(t,N) = F(t,n{t,x)), 

(7) 

which leads to 

ft — F t + F N n t , 

(8a) 


fx = FnU x . 

(8b) 

In particular, if we 

let f(t, x) = x and F(t, N ) = X(t, N ), then from Q we have 



0 = X t + X^riti 

1 XpjTl X ) 


which lead to 

x t 

q — n t — , 

Ajv 

(9a) 


-k = n x = —. 

A AT 

(9b) 
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Substituting (j9j) into (|8]), we obtain the transformation of derivatives in the Eulerian coordi¬ 
nates into those in the Lagrangian coordinates: 


ft 

fx 


. Xt 

hTT, 
A at 


’Jf 


X 


N 


(10a) 

(10b) 


We can repeatedly use (10) to transform higher-order derivatives from E-coordinate to 
L-coordinate. For example, if we denote g — ft and G = F t — then we have 


X t ^ _ X t _ . W ,,, _ XttXl- X?X* 


fu = g t = G t - = F tt - 2 F m XL + F NN (XXf - F N ^ tu N t ^ .(lOc) 

Aat Aat Aat Ajv 


Similarly, we have 


n _ ( rp t-i Xt X t ]\rX]\T — X t Xj\TN, 1 

Jtx — ln iA r — -CJVVTT--7^2- 

Aat A. n Aat 


fXX 


FnnXn — FnXnn 

x% 


(10d) 

(10e) 


In particular, when f(t,x ) = n(t,x ) and F(t,N ) = IV, we have q t = n tt = — Xtt 

n — —Jc, — n, — X tN X N -X t X NN i _iL _ _ X NN 

Hx l^t 'Hx ^3 5 dilU. Av# il'xx v-3 • 

. . . . N 

Similarly we can transform a variable and its derivatives from the Lagrangian coordinates 
to the Eulerian coordinates as follows: 


and 


F(t,N) 

F t 


F n 


f(t,X(t,N)), 


ft- fx— = ft + vf x , 
n x 


fx — 

n x 



(11a) 

(lib) 

(11c) 


F tt 

ftt ^vftx ”1“ ^ fxx “1“ i^t Wx) fxi 

(lid) 

FtN 

( ftx V fXX H“ '^xfx^T' 1 

k 

(lie) 

Fnn 

— (kfxx ^xfx)~j~^‘ 

(Ilf) 


In particular, when F(t, N) = X and f(t, x) = x, we have X t = v, Xn = — |, X tt = v t + vv x , 
X tN = -f, and X NN = -|f. 

In Table [lj where A = Vt and a = vt + vv x , we summarize transformations between L-P 
and E-S variables. In the table, variables are divided into three groups: the first group are 
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for k , q, v, the second group for simple derivatives of k and v in the Eulerian coordinates; and 
the third group for simple derivatives of X in the Lagrangian coordinates. The equivalent 
variables in E-P and L-S models can be easily obtained through changes of variables but are 
omitted here. 


Variables 

L-P 

E-S 

Density 

l 

Xn 

k 

Speed 

x t 

V 

Flow-rate 

Xt 

X N 

q 

Speed rate 

Xu - ^X tN 

v t 

Speed gradient 

± X tN 

Vx 

Speed curvature 

X-t.NN-X-N—X tN X NN 

X 3 

VXX 

Density rate 

XtNXff—Xt.X NN 

X 3 

k t 

Density gradient 

X-Nl V 

X N 

kx 

Acceleration rate 

X tt 

Vt + vv x 

Jerk 

x ttt 

Clt + V&x 

Speed difference 

x tN 

V x 

k 

Spacing difference 

Xnn 

kx 

k 3 


Table 1: Transformations between E-S and L-P variables 


2.2 Finite difference approximations of L-P derivatives 

In traffic flow, information travels in the increasing direction of N ; that is, the following 
vehicles have to adjust their behaviors according to the leading vehicles. Therefore, a simple 
upwind scheme can be used to numerically solve the L-P formulation as follows: first, we 
approximate the derivatives with respect to vehicle number by finite differences 


F N (t,N ) 


F{t,N) - F{t,N - AN) 
AN 


In particular, when AN = 1, a vehicle just follows its immediate leader, and we have 


F N (t,N) « F{t,N)-F(t,N- 1). (12) 

Similarly we can obtain the vehicle-discrete forms of speed and spacing differences 

X N (t,N) « X(t,N)-X(t,N-l), 

X tN (t, N) « X t (t,N)-X t (t,N-l), 

X NN (t,N ) « X(t,N) + X(t,N-2) -2X(t,N- 1 ). 
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Note that the right-hand sides of the approximations are discrete in vehicle number but 
still continuous in time. Such vehicle-difference L-P variables are used in time-continuous 
car-following models. 

Furthermore we can use the forward difference to replace the derivatives with respect to 
time in the vehicle-difference forms: 




F(t + At,N) - F(t, N) 
A t ' 


(13) 


where At should satisfy the CFL condition (Courant et al., 1928 Leclercq et al., 2007): 


At < max |0 , (S r )| with AN = 1, which can be derived from the L-S formulations. In particular, 
for the triangular fundamental diagram, we have At < Intuitively, information cannot 
pass for more than one vehicle during a time interval; i.e., a leading vehicle cannot impact 


the behavior of more than one vehicle, its follower, during a time interval. With both (12) 


and (13) we can obtain time-discrete car-following models from the L-P formulations. 


2.3 A two-way method 

Based on the transformations between E-S and L-P variables, we can then obtain equivalent E- 
S and L-P formulations when one of them is given. From the finite difference approximations 
of L-P derivatives, we can convert L-P formulations into car-following models and vice versa. 
Therefore, conversions between continuum and car-following models can be done in two steps, 
whose order depends on the direction of conversion. Note that the resulting E-S and L-P 
models are equivalent, but the car-following models are just approximately equivalent to the 
L-P models, subject to errors between derivatives and finite differences. 

Here we apply the method to derive equivalent time-continuous and time-discrete car¬ 
following models of the LWR model, the first-order continuum model, by discretizing (|5]): 

Xt(t,N) = e(X{t,N-l^X(t,N)), (14) 

X(t + At, N) = X(t,N) + Ate(X(t,N-l)-X(t,N)). (15) 


In particular, for the triangular fundamental diagram, we have 


X t (t,N) 
X(t + At, N ) 


min{rn, —( X(t,N — 1) — X(t,N ) — Sj)}, 

T 

X{t, N ) + At min-ftp, -(A(t, N — 1) — X(t, N) - S '■)}, 

T 


(16) 

(17) 


where r = —^- is the time gap, and Sj = 1/kj the jam spacing. Here the congested part of 
the time-continuous car-following model, (16), is essentially the same as the original Pipes 
car-following model, Equation 2.3 in (Pipes, 1953). If the minimum time At = r in the 


time-discrete Pipes model, (17), then we have 


X{t + r,N) = min{A(t, N) + tv f ,X(t,N — 1) — Sj}, (18) 
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which is Newell’s car-following model (Newell, 2002 Daganzo, 2006a). From the Hopf-Lax 
formula, we can show that the model is accurate subject to the assumption of uniform 
spacings between two consecutive vehicles (Evans, 1998). 


3 Continuum formulations of higher-order car-following 
models 

In this section, we apply the method in the preceding section to derive continuum formulations 
of higher-order time-continuous car-following models. In the first step, we approximate finite 
differences by partial derivatives with respect to vehicle numbers; in the second step, we 
convert L-P models to E-S continuum models. 


3.1 Second-order models 

In the literature, a large number of time-continuous car-following models have been developed 
to calculate a follower’s acceleration rate from its speed, spacing, and speed difference and 
can be written in the following form: 

X tt (t,N) = *(X t (t,N),AX(t,N),AX t (t,N)), (19) 

where vehicle N follows vehicle N — 1, the spacing A X(t, N) = X(t, N — 1) — X(t, N), and 
the speed difference A X t (t, N ) = X t (t, N — 1) — X t (t, N ). 

In the first step, from ( [l2| ), we can replace the discrete spacing AX(t, N ) by the continuous 
spacing — Xn, and the discrete speed difference AX t (t, N ) by the continuous speed difference 
—X t u- Then we obtain the corresponding L-P formulation of the car-following model, ( Jl9| ): 

Xtt = V{X t ,-X N ,-X tN ), (20) 

which is a second-order model as the highest order of derivatives is two. In the second step, 
from Table [l] we obtain the following E-S formulation: 

v t + vv x = $(i),^y), (21a) 

which, coupled with the conservation equation, 


h + (kv) x = 0, 

forms a second-order continuum model (PayneJ 1971 Whitham 1974). 
The following are some examples: 


(21b) 
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1. The linear General Motors model without delay, which was derived based on the 
stimulus-response principle (Pipes, 1953 Gazis et al., 1959), can be written as 


X a (t,N) = -A X,(t,N), 


( 22 ) 


where T > 0 is the reaction time. Thus T(n, |, ^) =_£jf, and the second-order 
continuum model is (we omit the conservation equation (21b) hereafter.) 


v t + ( v - ) v x = 0. 


(23) 


2. The nonlinear General Motors model without delay (Gazis et ah, 1961) can be written 
as: 


X u (t,N) = a 


X m AX t (t, N ) 

JaxmoF ’ 


where m and l are natural numbers, and a > 0. Thus T(u, = av m k 1 v x , and the 
second-order continuum model is 


v t + ( 


v — av 


m k l ~ 1 ) 


= 0 . 


(24) 


3. The Optimal Velocity Model (OVM) without delay |Bando et ah, 1995) can be written 


as: 


X tt (t,N) = T (0(AX(t,N))-X t (t,N)). 

Thus T(u, ^) = y(6*(r) — v), and the second-order continuum model is 


v t + vv x 


= 


(25) 


(26) 


Note that the model in (Newell, 1961), X t (t + T, N ) = 9(AX{t, N )), is equivalent to the 
Optimal Velocity Model if we approximate X t (t + T, N) by its Taylor series expansion: 
X t (t, N) + TX tt (t, N). 


4. The Generalized Force Model (GFM) (Helbing and Tilch 1998) can be written as: 

X a (t, N ) = i(0(AX(t, N)) - X t (t, N)) - N M Ax t(t, N )) e -[AX(t,N)-{d+TX t (uN))]/R\ 

where T is the acceleration time, T' < T the braking time, d the minimal vehicle distance, 
r the safe time headway (time gap), R' the range of the braking interaction, and ©(■) 
the ffeaviside function. Thus T(u, ^) = — v) — VxQ ^ k ^ e -[i/k-(d+Tv)]/R'^ an j 

the second-order continuum model is 




= 


( 27 ) 
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5. The Intelligent Driver Model (IDM) (Treiber et ah, 2000) can be written as: 

X tt (t, N) = a 


1 - 


X t (t,N) 


v f 


Xt(t,N)AX t (t,N) ' 


d + rX t (t,N) + ^^ 


A X(t, N) 


where a is the maximum acceleration rate, 6 the acceleration exponent, b the desired 
deceleration rate, r the safe time headway (time gap), and d the minimal vehicle 

distance. Thus ^(v, ^) = a 1 — — (^ + rkv + , and the second-order 

continuum model is 


v t + vv x = a 


v f 


1 - I — ) - ( dk + Tkv+ VV 


) 

‘l\fab ) 


(28) 


6. The Full Velocity Difference Model (FVDM) (Jiang et ah, 2001) is a combination of 


the linear General Motors model and the Optimal Velocity Model and can be written 
as: 


X tt (t,N) = -(d(AX(t,N))-X t (t,N)) + \AX t (t,N), 


(29) 


where T is the reaction time, and A the sensitivity coefficient. Thus T(n, ^) = 
b (0(4) — v) + A%, and the second-order continuum model is 


v t 




(30) 


Among the six continuum models, (23)-(30), the one for the linear General Motors model, 


(23), is a special case of the Aw-Rascle-Zhang model (Aw and Rascle, 2000 Zhang, 2002); 
the one for the Optimal Velocity model, (26), appeared in (Phillips, |1979 ); the one for the 
Full Velocity Difference Model, (30), is a special case of the Aw-Rascle-Greenberg model 


(Aw and Rascle, 2000| Greenberg, 2001); but, to the best of our knowledge, the other three 
models are new. 


3.2 Third-order models 

In some second-order car-following models, there is a delay term in the acceleration rate: 

X tt (t + r,N) = *(X t (t,N),AX(t,N),AX t (t,N)). (31) 

If approximating X tt (t + T', N ) by its Taylor series expansion, X tt (t, N ) + T'X ttt (t, N ), we 
then obtain a third-order car-following model: 

X m {t,N) = ^[*{X t {t,N),AX{t,N),AX t {t,N))-X tt {t,N)]. (32) 
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By using the same method in the preceding subsection, we can then derive the following 
third-order continuum model: 


h + (kv) x 
v t + vv x 

dt T va x 


0, 


a, 

1 

V 



1 

k’ 


Vx v 
k } 



(33a) 

(33b) 

(33c) 


where a is the acceleration rate, and a t + va x the jerk. 

For examples, we can derive such continuum models for the General Motors model with 


delay (Gazis et ah, 1959, 1961) and the Optimal Velocity Model with delay (Bando et al. 


1998). Further higher-order models can also be derived by using this method. 


4 Car-following formulations of higher-order continuum 
models 

In this section, we apply the conversion method in Section 2 to derive car-following models 
from higher-order continuum models. 


4.1 A review of higher-order continuum models 

We give a brief review of existing second-order continuum models and categorized them 
into four families. All second-order models consist of the conservation equation, (21b), and 
an additional equation in speed. Here we just list the equation in speed and denote the 
speed-density relation by r)(k) = 9(1/k): 


1. The Phillips family: 

• ([Phillips, 1979): v t + vv x = ^(r)(k) — v)] 

• (Ross, 1988): v t + vv x = ^(vf — v ), which was shown to be unrealistic in (Newell 


1989); 


Liu et al. 
unction of 


1998): Vt = Tjkjivik) — v), in which the relaxation time T(k) is a 
density k. 


2. The Payne family: 


(Payne, 

1971) 

: v t + w x 2 Tkk x =T(ri(k) v); 

(Whitham, 

1974 

): v t + vv x + ^-k x = ^( rj(k ) - u); 

(Phillips 

1971 

)■ 

v t + vv x + ^p-k x = if(r](k) — v ), where p(k) is the traffic pressure 


function; 
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• dKuhne , |1987 j Kuhne and Michalopoulos 1992; Kerner and Konhauser 1993): 
v t + vv x + ^fk x = rj(k ) — v) + f -f-v xx , which has a viscous term; 

• (Zhang, 1998): v t + vv x + {k ' n [ k)) ' k x = £( rj(k ) - v). 

3. The Aw-Rascle family: 


(Aw and Rasclc, 2000): 


v t + (■v - kp'(k))v x = 0; 


(34) 


(Zhang, 2002): 


v t + (v + krf{k))v x = 0; 


(35) 


(Greenberg, 

2001): 

(Jiang et al. 

2002) 


v t + (v- c 0 )v x = —{q{k) - v)- 


(36) 


• (Xue and Dai, 2003): v t + (v + k-^^r]'{k))v x = (jj{k) — v ), where t r is driver’s 
reaction time. 

4. The Zhang family: 

• (Zhang, 2003): v t + (v + 2 f3c(k))v x + ^ —k x = rj(k ) — v) + n(k)v xx , where 

c(k) = krj’Jk), n(k) = 2/3Tc 2 (k). 

From these models we have the following observations: 

• Models of the Phillips family are the simplest: those of the Payne-Whitham family 
have one extra term with the density gradient k x , those of the Aw-Rascle family have 
one extra term with the speed gradient v x ; and those of the Zhang family have both. 

• Most models have the total derivative in v, v t + vv x , and some just have the partial 
derivative, v t . 

• Some models have the relaxation term, j.(j}(k) — v ), and some do not. 

• Some models have the viscous term v XXl but most do not. 

In the literature, third- or even higher-order continuum models have been proposed 


(Phillips, 1979 Helbing, 1995, 1996 Colombo, 2002a b). But in this study we focus on the 
second-order models. 
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4.2 Car-following models of second-order continuum models 

Given a continuum model, in the first step, we convert them into L-P formulations by 
following Table [I] in the second step, we approximate partial derivatives with respect to 
vehicle numbers by finite differences. In this subsection we only use the car-following models 
of the Aw-Rasclc family as examples. 

Even though the continuum formulations of General Motors and Full Velocity Difference 
models belong to the Aw-Rascle family, we can obtain more general car-following models of 
the following general Aw-Rascle model: 


v t + {v- kp'(k))v x = (v(k) - v). 

First, from Table [I] we can have the L-P formulation of (37): 

1 


X tt = 


t(—i/x n ) 


(d(-X N )-X t )-p'(-l/X N ) 


X 


tN 


X 2 n 


(37) 


(38) 


Further replacing —X N by A X(t,N) and —X tN by A X t (t,N), we obtain the following 
car-following model: 


X tt {t,N) = 


T(l/AX(t,N)) 


mX(t, N)) - X,(t, N)) + p\l/AX{t, N)) X^X£ 39) 


The new car-following model (39) combines a generalized Optimal Velocity model (the first 
term on the right-hand side) and a nonlinear General Motors model (the second term on the 
right-hand side). This model is also in the form of (19). 


For example, the car-following formulation of the Aw-Rascle-Zhang model, (35), is 

A X t (t,N) 


X tt (t, N) = -rf(l/AX(t,N)) 


AX(t, A) 2 ’ 


(40) 


which is different from X u (t, N ) = • The latter was used in (Zhang 

derive (35) by Taylor series expansions. As another example, the car-following, (36 , has the 
following car-following formulation: 

X u (t,N) = l(g(A X(t, TV)) - X t (t, N)) + , (41) 

which is a combination of the optimal velocity model and a nonlinear General Motors model. 
Note that (36) was derived with Taylor series expansion from the FVD model, which is a 
combination of the optimal velocity model and the linear General Motors model. 

To the best of our knowledge, both (40) and (41) and their corresponding continuous L-P 
formulations are new. Actually from (40), we can choose different pressure functions, p(k), 
to obtain other new car-following models. Following the same procedure, we can obtain more 
new car-following models from continuum models of other families. 


2002) to 
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5 Equivalence between continuum and car-following 
models 

Different from the first-order models, whose analytical properties and numerical errors have 
been well understood, higher-order models are much more challenging to analyze, as there 
exist no systematic methods for solving higher-order E-S or L-P models. Thus instead of 
determining the conditions when continuum and car-following models are equivalent, in 
section we attempt to examine their equivalence with respect to steady-state speed-density 
relation and stability of steady states. 


5.1 Steady-state speed-density relations 


In this study we adopt the following definitions of steady states: for car-following models, 
Xt(t,N) = v and A X(t,N) = s = i.e., all vehicles have the same constant speed v and 
spacing s; for continuum models, k(t,x ) = k, and v(t,x) = v; i.e., the density and speed are 
constant with respect to both time and location. 


For a general second-order car-following model in (19), we have X t t(t,N ) = 0, and 


A X t (t, N) = 0 in steady states. Thus by solving the following equation: 


0 ) = 0 , 


(42) 


we can obtain the speed-density relation, v = rj(k). For the corresponding continuum model 


in (21), we have v t = 0, and v x = 0. Thus the speed-density relation also satisfied (42). 


That is, the car-following and continuum model share the same speed-density relation and, 
therefore, the same fundamental diagram. 

In the following we examine three types of higher-order models: 

1. Most of the car-following models in Section 3 and most of the continuum models 
reviewed in Section 4.1 have a simple relaxation term, 4(0(AA"(f, N)) — X t (t,N )) or 
4 [rj[k) — v), and their steady-state speed-density relations are simply v = r}(k) = 9(i)- 


2. For the IDM and its continuum counterpart, (28), no simple relaxation term is given. 


We have from (42) 


v 


— — (dk + rkv) = 0, 

v f, 


which leads to the following density-speed relation, k = ij 1 (n): 

V2 

d + TV \ \Vf ) 


k = 
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Thus we have the following observations: (i) v = Vf if and only if k — 0; (ii) k decreases 
in v, and 0 < k < A (iii) 0 < v < Vf when k > 0. In particular, when S —>■ oo, we have 


v 


min 




which corresponds to the triangular fundamental diagram. The speed-density relation 
is consistent with that derived from the car-following model in (Treiber et ah, 2000). 


3. For the original Aw-Rascle model, (34), the corresponding car-following model is given 
by 


X tt (t,N) = p\l/AX(t,N)) 


A X t (t,N) 
AX(t, N) 2 ' 


for which the steady states occur if and only if X t {t , N — 1) = X t (t, N ), regardless of 
the spacing. Therefore, the speed-density relation in steady states is not unique. This 
problem also appears in the General Motors car-folio wing model ( Nelson[ 1995). It was 


argued in (Gazis et ah, 1959) that the steady-state solutions fall on a valid, unique 


fundamental diagram if the initial conditions satisfy the same fundamental diagram. 
However, in a road network where vehicles constantly join and leave a traffic stream at 
merging and diverging junctions as well as on a multi-lane road, the steady-state speed- 
density relation cannot be maintained in such models. A simple correction to these 
models is to add a relaxation term in both continuum and car-following formulations. 


Motors model, (22) 


For example, the FVD model, (29), can be considered a correction to the linear General 


5.2 Stability of steady states 

In this subsection we analyze the stability of the general second-order car-following model, 
(19), and the corresponding continuum model, (21). 

For (19), we follow (Chandler et ah, 1958) to analyze the string stability of a platoon of 
vehicles and have the following result. 


Theorem 5.1 The car-following model, (19), is string stable if and only if 

*S(t«o, So, 0) > 21,(»o, so, o). 


(43) 


Note that the string stability condition, (43), is irrelevant to the partial derivative with respect 
to the speed difference. 


Proof. Assuming that the leading vehicle’s speed X t (t,N — 1) = no + cn- ie tu>t and the 
following vehicle’s speed X t (t, N ) = uo + ejve ia ' t are small monochromatic oscillations around a 
steady state with a speed v 0 and spacing s 0 , we then have X u (t, N ) = iuje^e lwt , X(t,N — 1) = 
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v 0 t + e -^e iu)t + X(0, N - 1), and X(t, N) = v 0 t + e ^e iut + X(0, N ). Thus, A X(t, N) = 
s 0 + £N -}~ ejv e lut , where s 0 = -V(0, N — 1) — A(0, N ) is the initial spacing, and A X t (t, N ) = 


(cat-i — 6N)e lU3t . From (19), we have 

tue N e iut = y(vo + e N e- t ,s 0 + eN - 1 - eN ^ t 


IU) 


e , (ov-i ■” 


With the Taylor series expansion of the right-hand side, we then have 


lue^e 


,iu)t 


iut i YTr r, n\ 1 C-N iujt 


^(Vo, s 0; 0) + T„(fo, So, + ^(uo, So, 0) 


ICO 


e lU1 + s 0) 0)( e iV-l 


Since v 0 and s 0 satisfy T(n 0 , Sq, 0) = 0 in steady states, we have 


IUJ6n 


^v{vq,s 0 , 0)e N + ty s (vo,s o ,0 )-^.—— + l / Ai)(ro, s 0 , 0)(ejv-i - cat)- 


IU 


Thus we have the ratio of over ejv-i as 

e N ^ s (v o ,s o ,0) +iujty Av (v o ,s o ,0) 


Cat-1 


-uj 2 - iu^ v (v 0 , So, 0) + T s (u 0 , So, 0) + iuj^ Av (v 0 , s 0 , 0)' 


Then the car-following model, (19), is string stable if and only if | ctv/ctv—i | < 1 for any 
frequency cu; i.e., the amplitude of oscillations decays along a platoon. Hence the string 
stability condition can be written as So, 0) + co 2 > 2T s (n 0 ,s 0 ,0) for any uj, which is 

equivalent to (43). ■ 

For example, for the Optimal Velocity model, (25), the string stability condition, (43), 
can be simplified as 


2T9 s (s 0 ) < 1. 


In particular, for the triangular fundamental diagram, 0(s o ) = min{i,y, y(so — Sj)}, and we 
have the following conclusion: (i) it is always stable in uncongested traffic, since 9 s (s 0 ) = 0; 
(ii) in congested traffic, 9 s (s 0 ) = y, and the stability condition is equivalent to 


i.e. 


the relaxation time is less than half of the time-gap. 


For (21), we follow (Whitham 1974, Chapter 3) to analyze its linear stability. We first 
linearize it about a steady state at k 0 and v 0 , which satisfy T(n 0 , yy,0) = 0, by assuming 
k = k 0 + e and v = v 0 + £, where e(x,t) and £(x,t) are small perturbations. Omitting 
higher-order terms, we have (21) 


O + HjO + ko9, x — 0, 

ft + (u 0 - ^Ai;(uo,so,0)so)^ = ^v(vo,s 0 ,0)f - tf a (u 0 ,so,0)sjje, 


- e N )e^. 
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where Sq = Taking partial derivatives of the second equation with respect to x and 
substituting from the first equation we obtain 

d d \ f d d\ 

— + (v 0 - ^Av(vo, So, °) s °)^: ) ( -Q f + v 0 fa ) e = ^v( v o, So, 0)(e t + v 0 e x ) + T s (u 0 , s 0 , 0)s 0 e*. 

With the exponential solution e = where m is a real number, then the linearized 

system is stable if and only if two solutions of uj in the following equation have negative 
imaginary parts: 

[cu - (u 0 - ^Av(v 0 , s 0 , 0)s 0 )m ](u - v 0 m) + s 0 ,0)(-o; + v 0 m ) + T s (u 0 , s 0 ,0)s 0 m] = 0, 

which can be re-written as co 2 + (2b\ + i26 2 )w + (di + id 2 ) =0 with 

26i = — (2u 0 - T A t,(u 0 , s 0 , 0)s 0 )m, 2b 2 = -^ V (v 0 , s 0 ,0), 
di = (v 0 - * Av (v 0 , s 0 , 0)s 0 )u 0 m 2 , d 2 = (dh(u 0 , s 0 , 0)u 0 + ^ s (^o, s 0 ,0)s o )m. 


According to (Abeyaratne 2014), both roots of the equation have negative imaginary parts if 


and only if b 2 > 0 and 4 bib 2 d 2 — Adib\ > d\, i.e., \& v (i? 0 , s 0 ,0) < 0, and 
T 2 (u 0 , s 0 , 0) + T„(u 0 , s 0 , 0)T s (u 0 , s 0 , 0)^ Av (v o , s 0 , 0) < 0. 

Thus for the optimal velocity model, the linear stability condition becomes 4/ 2 (uq, s 0 , 0) < 0, 


which is impossible. Such an inconsistency was observed in (Abeyaratne, 2014), where a 


modified version of (21). However, the modification may not apply to general models. Such 


an inconsistency could be interpreted in two ways: first, it is possible that string stability is 
different from linear stability; second, it is possible that the continuum and car-following 
models obtained with the conversion method may not be equivalent, due to the upwind 


difference method for discretization of the L-P models, (12). The real underlying reason is 


subject to future exploration and seem to be quite challenging. 


6 Conclusion 


Inspired by recent studies on equivalent formulations of the LWR model (Daganzo 2006a|b 


Leclercq et ah, 2007; Laval and Leclercq, 2013), in this study we presented a unified 


approach to convert higher-order car-following models into continuum models and vice 
versa. The conversion method consists of two steps: equivalent transformations between 
the secondary Eulerian (E-S) formulations and the primary Lagrangian (L-P) formulations, 
and approximations of L-P derivatives with finite differences. In a sense, continuum and 
car-following models are equivalent, subject to the errors between finite differences and 
derivatives. The conversion method is different from the traditional method based on Taylor 
series expansion, as it works in both directions. We used the method to derive continuum 
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models from general second- and third-order car-following models and derive car-following 
models from second-order continuum models. We demonstrated that higher-order continuum 
and car-following models have the same fundamental diagrams. We also derived the string 
stability conditions for general second-order car-following models and the linear stability 
conditions for general continuum models. But we found that that they are different. The 
underlying reason of such inconsistency is subject to future research, and the conditions when 
continuum and car-following models are equivalent are yet to be determined. This study 
further highlights substantial challenges associated with higher-order models. 

Through this study, we revealed underlying relationships among many existing models. 
For example, the Optimal Velocity car-following model (Bando et ah. 1995) and the Phillips 
second-order continuum model (Phillips, 1979) are equivalent with the proposed conversion 
method. In addition, with the conversion method we derived a number of new continuum 
models from existing car-following models and new car-following models from existing 
continuum models. Furthermore, as for the LWR model in (Laval and Leclcrcq, 2013 ), we 
can extend the method to derive six equivalent continuous formulations of a traffic flow model: 
the primary and secondary formulations in (t,x)~, (' t,N )-, and (x, N) -coordinates. Different 
formulations can be helpful for solving different initial-boundary value problems. However, 
we have yet to develop simple and systematic tools to analyze and solve such higher-order 
models, as traditional tools for Hamilton-Jacobi equations and hyperbolic conservation laws 
cannot be directly applied. 

In the future, it is possible to extend the conversion method for more complicated traffic 
flow models, for example, with multiclass drivers and stochastic driving characteristics, 
with lane-changing, merging, and diverging behaviors. We can develop hybrid models from 
different formulations. We can use different data sources to calibrate a model and predict 
traffic conditions. Finally, we can develop traffic control and management strategies with 
different formulations. 
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